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We present measurements of differential cross sections in inclusive Z/-y* plus jet production in a data sample 
of 1 fb^ 1 collected with the D0 detector in proton antiproton collisions at y/s = 1.96 TeV. Measured vari- 
ables include the Z/7* transverse momentum (p z ) and rapidity (y z ), the leading jet transverse momentum 
(p^ et ) and rapidity (j/-' et ), as well as various angles of the Z+jet system. We compare the results to different 
Monte Carlo event generators and next-to-leading order perturbative QCD (NLO pQCD) predictions, with 
non-perturbative corrections applied. 



X 
CD 

Oh' 

CD 
43 



(N 
> 
O 
oo 
m 



On 
O 



X 



1. Introduction 

To make discoveries at the Tevatron and the LHC, 
background processes will need to be measured and 
simulated with a level of accuracy that will be com- 
parable to the signficance of those new physics sig- 
nals. There are several programs on the market that 
can simulate hadronic interactions at next-to-leading 
order (NLO) accuracy, but the processes included in 
these programs are limited. Matrix element plus par- 
ton shower (MEPS) programs simulate a more com- 
prehensive set of processes, typically at leading-log 
(LL) or leading order (LO), and rely on models to 
simulate emissions and fragmentation associated with 
higher order processes. These programs have been 
employed regularly for background simulation at the 
Tevatron in recent years, notably in the Higgs searches 
[l[ and the discovery of the production of single top 
quarks @. 

Measurements of Z/j* + jets processes are valuable 
for two principle reasons. Z/"f* production provides 
a hard scale, which, along with associated jet produc- 
tion, is an ideal environment to test perturbative QCD 
(pQCD). The leptonic decay of the Z/-f* provides 
a clean signal for reconstruction of the events, and 
small background contamination. The test of pQCD is 
made by comparing the measurements to NLO pQCD 
predictions. Z /~/* + jets also makes up a major back- 
ground of many new physics searches at both the 
Tevatron and LHC. Therefore, these data measure- 
ments unfolded to the particle level are useful for tun- 
ing LO simulation programs which are heavily relied 
upon to model background processes. 

The Tevatron measurements presented here of 
Z/l* + jets differential cross sections are compared to 
predictions by NLO pQCD in MCFM i, MEPS pro- 
grams ALPGEN (1 and SHERPA jf, and PS pro- 
grams HERWIG @ and PYTHIA_[7L_ The measure- 
ments have either been published [g| ^| or have been 
submitted for publication [l(| at the time these pro- 
ceedings were written. ALPGEN employs the MLM 
algorithm to ensure jets originating from the ma- 
trix element and the parton shower are not double 



counted. SHERPA is a CKKW-inspired model which 
uses a reweighting of the matrix elements to achieve 
the same appropriate jet configurations. A detailed 
description of these programs can be found in [llj . 



2. Data Selection 

The measurements are made with the D0 detector, 
which is described in detail elsewhere [13]. The 
analyses were performed in the Z/j* — + fifi and 
Z/j* — > ee decay channels. Z/j* — > fifi focused on 
the Z+l jet inclusive events, and differential cross 
sections were made in a variety of variables, includ- 
ing many angular variables involving the decay ob- 
jects. The Z/j* — > ee analysis focused on differen- 
tial cross section measurements as a function of p^ e ' , 
in the 1, 2 and 3 jet inclusive samples. In both 
channels, the data are corrected to particle level to 
eliminate the effects arising from detector resolution 
and efficiency. In order to keep these corrections 
small, a limited particle level phase space for the 
measurements was chosen that corresponds closely to 
the detector level selection cuts. Some effort was 
made to keep the particle level phase space in the 
Z/j* — > fifi and Z/j* — > ee channels close, but some 
differences exist due to the different nature and detec- 
tor manifestations of electrons and muons. At particle 
level, the Z/j* is reconstructed using the two highest 
Pt leptons, and the jets are reconstructed using the 
D0 Runll Midpoint Cone algorithm with a cone ra- 
dius R=0.5. In the Z/j* — ► fj.fi analysis, the particle 
level phase space is defined with the following restric- 
tions: 



• 65 < < 115 GeV 

• p 3 T et > 20 GeV 

• \y jet \ < 2.8 

• |y M | < 1.7 

In the Z/j* — > ee analysis, the particle level phase 
space is defined with the following restrictions: 
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• 65 < M ee < H5 GeV 

• P J T et > 20 GeV 

• \y jet \ < 2.5 

3. Theoretical Predictions 

Due to the evolving nature of theoretical predic- 
tions, the programs used to compare to the data are 
not consistent in all distributions. This is due to the 
fact that the data are published at different times, 
and effort is made to always compare to the most up- 
to-date theoretical predictions. Next-to-leading order 
perturbative QCD (NLO pQCD) predictions are cal- 
culated using the program MCFM; parton-to-hadron 
corrections are applied to the MCFM predictions, as 
calculated by PYTHIA. 

The meausurements of Z/^* — ► \x\x + jets are com- 
pared to different predictions than those compared in 
the Z/"f* — > ee + jets measurements. The programs 
and their version numbers are summarized in Table UJ 



theoretical uncertainties. SHERPA generally gives 
the best and most consistent description of the an- 
gular distributions, although the normalization of the 
data is not reproduced. The description of the data 
by the LO MC models is mixed. All LO MC mod- 
els suffer from large scale uncertainties which impair 
their ability to make precise predictions. 

5. Conclusions 

Several differential cross sections of Z/j* + jet + X 
events measured with the D0 detector have been pre- 
sented. The data are generally consistent with predic- 
tions from NLO pQCD, although some LO programs 
can reproduce the shape of the data better than NLO, 
due either to their inclusion of higher parton multiplic- 
ity matrix elements than can be currently included in 
a fixed order pQCD calculation, or an optimized tune 
of the Monte Carlo. These data should be useful for 
continued tuning of these and other Monte Carlo pro- 
grams. 



Table I Theoretical predictions and their versions for 
Z/-y* — > and Z/'y* — * ee + jets data comparisons. 



Program 


Z/y* -> nn 


Z/j* — > ee 


MCFM 


5.4 


5.3 


ALPGEN 


2.13 


2.13 


SHERPA 


1.1.3 


1.1.1 


PYTHIA 


6.420 


6.416° 


HERWIG 


6.510 


6.510 



"PYTHIA 6.325 is used for interfacing ALPGEN with parton 
showers 



4. Results 

Particle level differential cross sections as func- 
tions of a variety of variables are shown in Fig- 
ures 1-13 for Z/j* + jets events. Figures 1-4 
and 8-13 illustrate measurements that were made 
in the Z '/j* — > decay channel, while Figures 
5-7 show measurements that were made in the 
Z/'y* — > ee decay channel. Figures 8, 10 and 12 
consider Z/-f* + jets events in the restricted range 
Z/7* > 25 GeV, while Figures 9, 11 and 13 re- 
strict this range further to Z/j* > 45 GeV. These 
restrictions are made to reduce the dependence on 
the underlying event and multiple interactions. All 
if/7* — > nn results are made with Z/j* + 1 jet in- 
clusive samples. The Z/j* — > ee results examine the 
Z/'y* + 1, 2 and 3 jet inclusive samples. A good de- 
scription of the data by NLO pQCD is found in all dis- 
tributions, taking into account the experimental and 
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Figure 1: The measured cross section in bins of leading 
for Z/j* + jet + X events, and predictions from NLO 
pQCD and SHERPA are shown in the upper left plot. 
The ratio of data to SHERPA are shown in the lower left 
plot. The ratio of LO and NLO MCFM predictions to 
SHERPA and associated scale and PDF uncertainties are 
also shown in the lower left plot. The ratio of PYTHIA 
and HERWIG+JIMMY to SHERPA are shown in the 
upper right plot. The ratio of ALPGEN+PYTHIA and 
ALPGEN+HERWIG to SHERPA are shown in the lower 
right plot. The ratio of data to SHERPA are included in 
all plots for comparison. 
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Figure 2: The measured cross section in bins of leading 
y z for Z/y* + jet + X events, and predictions from NLO 
pQCD and SHERPA are shown in the upper left plot. 
The ratio of data to SHERPA are shown in the lower left 
plot. The ratio of LO and NLO MCFM predictions to 
SHERPA and associated scale and PDF uncertainties are 
also shown in the lower left plot. The ratio of PYTHIA 
and HERWIG+JIMMY to SHERPA are shown in the 
upper right plot. The ratio of ALPGEN+PYTHIA and 
ALPGEN+HERWIG to SHERPA are shown in the lower 
right plot. The ratio of data to SHERPA are included in 
all plots for comparison. 
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Figure 3: The measured cross section in bins of leading 
p^ e * for Z/-y* + jet + X events, and predictions from NLO 
pQCD and SHERPA are shown in the upper left plot. 
The ratio of data to SHERPA are shown in the lower left 
plot. The ratio of LO and NLO MCFM predictions to 
SHERPA and associated scale and PDF uncertainties are 
also shown in the lower left plot. The ratio of PYTHIA 
and HERWIG+JIMMY to SHERPA are shown in the 
upper right plot. The ratio of ALPGEN+PYTHIA and 
ALPGEN+HERWIG to SHERPA are shown in the lower 
right plot. The ratio of data to SHERPA are included in 
all plots for comparison. 




Figure 4: The measured cross section in bins of leading 
y jet for Z/j* + jet + X events, and predictions from NLO 
pQCD and SHERPA are shown in the upper left plot. 
The ratio of data to SHERPA are shown in the lower left 
plot. The ratio of LO and NLO MCFM predictions to 
SHERPA and associated scale and PDF uncertainties are 
also shown in the lower left plot. The ratio of PYTHIA 
and HERWIG+JIMMY to SHERPA are shown in the 
upper right plot. The ratio of ALPGEN+PYTHIA and 
ALPGEN+HERWIG to SHERPA are shown in the lower 
right plot. The ratio of data to SHERPA are included in 
all plots for comparison. 



Proceedings of the DPF-2009 Conference, Detroit, MI, July 27-31, 2009 



■5 



7l-i (-»ae) + 1j8t + X 
65<M aB <115GeV 

R™ = 0.5, | <2.5 




20 30 40 50 



)0 200 300 

P T (1 st jet) [GeV] 



-♦-Data -■-PYTHIASO 

— HERWIG+JIMMY - -Scale unc. 

— PYTHIAQW 
(c) — Scale unc. 




-t-Data SHERPA 
— ALPGEN+PYTHIA 
2.0 - — Scale uric. 




20 30 40 50 



10 200 300 

p T (1 s, jet) [GeV] 



y 3.d 

3 2 C 



DO Run II, L=1 .04 fb 1 -*- Data at particle level 
^ MCFM LO 

- (a) S *w 



Z/y [-* ee] + 3 jets + X 
65<M BO <115GeV 

R™ = 0^5, |y' e '|<2.5 



- Data 

- MCI M _0 



40 50 60 

p T (3"|et) [GeV] 



-Data -.-PYTHIASO 
-HERWIG+JIMMY Scale unc. 

PYTHIA QW 




p,(3 |et) [GeV] 



Figure 5: The measured cross section for Z/j* + 1 jet + 
X events in bins of leading p 3 ^ 1 . Predictions from NLO 
pQCD are compared to the data in the upper left plot. The 
ratio of data and several LO programs to MCFM NLO are 
shown in the other plots. MCFM NLO predictions agree 
with the data. The data can be described by a selection 
of the LO programs, although there is a lot of freedom in 
the predictions due to the choice of PDF, renormalization 
and factorization scale, tune and underlying event model. 



Figure 7: The measured cross section for Z/j* +3 jet + 
X events in bins of the third leading p^f . Predictions 
from NLO pQCD are compared to the data in the upper 
left plot. The ratio of data and several LO programs to 
MCFM NLO are shown in the other plots. MCFM NLO 
predictions agree with the data. The data can be described 
by a selection of the LO programs, although there is a lot 
of freedom in the predictions due to the choice of PDF, 
renormalization and factorization scale, tune and underly- 
ing event model. 




Figure 6: The measured cross section for Z/j* + 2 jet + 
X events in bins of the second leading p 1 ^ 1 . Predictions 
from NLO pQCD are compared to the data in the upper 
left plot. The ratio of data and several LO programs to 
MCFM NLO are shown in the other plots. MCFM NLO 
predictions agree with the data. The data can be described 
by a selection of the LO programs, although there is a lot 
of freedom in the predictions due to the choice of PDF, 
renormalization and factorization scale, tune and underly- 
ing event model. 
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Figure 8: The measured cross section in bins of A$ be- 
tween the Z/'y* and leading jet for Z/y* + jet + X events 
with pf larger than 25 GeV. Predictions from NLO pQCD 
and SHERPA are compared to the data in the upper 
left plot. The ratio of data and predictions from NLO 
pQCD, ALPGEN, HERWIG and PYTHIA to the pre- 
diction from SHERPA are shown in the lower plots. The 
shape of the distribution is best described by SHERPA, 
but there is a normalization disparity. 



Figure 10: The measured cross section in bins of Ay be- 
tween the Z/~/* and leading jet for Z/y* + jet + X events 
with pf larger than 25 GeV. Predictions from NLO pQCD 
and SHERPA are compared to the data in the upper 
left plot. The ratio of data and predictions from NLO 
pQCD, ALPGEN, HERWIG and PYTHIA to the pre- 
diction from SHERPA are shown in the lower plots. 
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Figure 9: The measured cross section in bins of A$ be- 
tween the Z/y* and leading jet for Z/"/* + jet + X events 
with pf larger than 45 GeV. The restricted pf range 
is chosen to isolate the measurement from biases due to 
the underlying event. Predictions from NLO pQCD and 
SHERPA are compared to the data in the upper left plot. 
The ratio of data and predictions from NLO pQCD, ALP- 
GEN, HERWIG and PYTHIA to the prediction from 
SHERPA are shown in the lower plots. 
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Figure 11: The measured cross section in bins of Ay be- 
tween the Z/~/* and leading jet for Z/y* + jet + X events 
with pf larger than 45 GeV. The restricted pf^ range 
is chosen to isolate the measurement from biases due to 
the underlying event. Predictions from NLO pQCD and 
SHERPA are compared to the data in the upper left plot. 
The ratio of data and predictions from NLO pQCD, ALP- 
GEN, HERWIG and PYTHIA to the prediction from 
SHERPA are shown in the lower plots. 
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Figure 12: The measured cross section in bins of yboost 
for Z/-y* + jet + X events with pf larger than 25 GeV, 
where yboost is the average rapidity of the reconstructed 
Z boson and the leading jet rapidity. Predictions from 
NLO pQCD and SHERPA are compared to the data in 
the upper left plot. The ratio of data and predictions 
from NLO pQCD, ALPGEN, HERWIG and PYTHIA 
to the prediction from SHERPA are shown in the lower 
plots. The shape of the distribution is best described by 
SHERPA, but there is a normalization disparity. 



Figure 13: The measured cross section in bins of yboost 
for Z/'y* + jet + X events with larger than 45 GeV. 
Predictions from NLO pQCD and SHERPA are compared 
to the data in the upper left plot. The ratio of data and 
predictions from NLO pQCD, ALPGEN, HERWIG and 
PYTHIA to the prediction from SHERPA are shown in 
the lower plots. 



